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Differential rotation of r-modes is investigated within the nonlinear theory up to second order in 
the mode amplitude in the case of a slowly-rotating, Newtonian, barotropic, perfect-fluid star. We 
find a nonlinear extension of the linear r-mode, which represents differential rotation that produces 
large scale drifts of fluid elements along stellar latitudes. This solution includes a piece induced 
by first-order quantities and another one which is a pure second-order effect. Since the latter is 
stratified on cylinders, it cannot cancel differential rotation induced by first-order quantities, which 
is not stratified on cylinders. It is shown that, unlikely the situation in the linearized theory, r- 
modes do not preserve vorticity of fluid elements at second-order. It is also shown that the physical 
angular momentum and energy of the perturbation are, in general, different from the corresponding 
canonical quantities. 
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I. INTRODUCTION 

R-modes in Newtonian gravity were first studied more than twenty years ago P, 0, 0j- In recent years, it was 
discovered numerically 0, and soon afterwards confirmed analytically that these modes, analogous to Rossby 
waves in the Earth's atmosphere and oceans, are driven unstable by gravitational radiation reaction in perfect fluid 
stars with arbitrary small rotation. 

The interest in the astrophysical implications of r-modes increased dramatically when it was shown [|| Q that in 
a newly born, hot, rapidly-rotating neutron star the radiation reaction force dominates bulk and shear viscosity for 
enough time to allow most of the star's angular momentum to be radiated away as gravitational waves. As a result, 
the neutron star spins down to just a small fraction of its initial angular velocity, thus providing a possible explanation 
for the relatively small spin rates of young pulsars in supernova remnants. For typical equations of state of a neutron 
star it was estimated 0,0 that the r-mode instability spins down a young neutron star to a period of about 10-20 ms. 
This is comparable to the inferred initial periods Pq of the fastest pulsars associated with supernova remnants, such 
as the Crab pulsar PSR B0531+21 (P = 19 ms), or PSR J0537-6910 (P < 16 ms for braking index n = 3 8]). The 
detectability of gravitational waves from such sources have been analyzed within a simple phenomenological model for 
the evolution of the r-mode instability, with the conclusion that gravitational waves emitted from an young, rapidly 
rotating neutron star in the Virgo cluster could be detected by enhanced versions of laser interferometer detectors 
0- However, recent results on the nonlinear saturation of the r-mode energy indicate that enhanced LIGO detectors 
could detect gravitational waves from these sources only up to a distance of about 200 kpc ^(j • 

Another interesting astrophysical implication of r-modes is related to the fact that gravitational- wave emission due 
to mode instability could balance the spin-up torque due to accretion of neutron stars in low-mass X-ray binaries 
|lll Il2| , thus limiting the maximum angular velocity of these stars to values consistent with observations (recent 
results indicate that the fast end of the spin distribution could be around 600 Hz, well below the theoretical maximum 
[l3j|). If the viscosity of the neutron star is a decreasing function of temperature, then the star undergoes a cyclical 
evolution of spin-up due to accretion and spin-down due to gravitational radiation emission [T^j . Assuming that the 
saturation amplitude of the r-mode is of order unity, it was estimated that the time spent by the star in the 
spin-down phase is about one year, just a small fraction of the full duration of the cycle (several million years), thus 
making it unlikely that any such source in our Galaxy is presently emitting gravitational waves. However, recent 
results indicate that the saturation amplitude of the r-mode could be much smaller than unity, implying that the 
duration of the spin-down phase could be a much larger fraction of the cycle |l5j | ; in that case, gravitational radiation 
from such sources in our galaxy could be detected by enhanced LIGO detectors |ic| . 

While initial research on r-modes was carried out within the linearized theory, a deeper understanding of these 
modes and its astrophysical implications requires considering the nonlinear theory. 
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An important nonlinear effect, that has been investigated recently by several authors, is differential rotation induced 
by r-modes. This differential rotation could wind up a magnetic field in an old, accreting neutron star in an X-ray 
binary, leading to a gamma-ray burst . Differential rotation could also interact with the magnetic field of a newly 
born, rapidly rotating neutron star, limiting the growth of the r-mode instability or, for strong magnetic fields, even 
preventing it from developing [171 Il8l Il9| . 

An analytical expression for differential rotation induced by r-modes was first derived using the linearized fluid 
equations by expanding the velocity of a fluid element located at a certain point in powers of the mode amplitude, 
averaging over a gyration, and retaining only the lowest-order nonvanishing term |l7| . Since this procedure is not 
equivalent to solving the nonlinear fluid equations, the analytical expression thus obtained is just an approximate one. 
Differential rotation was also reported in a useful toy model of a thin spherical shell of a rotating incompressible fluid 
|20|. However, in this case, differential rotation is of a different nature than the one referred to above; it is driven by 
radiation reaction. The existence of differential rotation related to r-modes was confirmed by numerical simulations 

In this paper, differential rotation induced by r-modes is investigated within the nonlinear theory up to second order 
in the mode amplitude a in the case of a Newtonian, barotropic, perfect-fluid star rotating with constant angular 
velocity fl. Second-order quantities are obtained by expanding simultaneously in powers of a and Q. In order to 
consistently neglect higher-order terms arising from the expansion in ft, it is assumed that a S> (ft /flu) 2 , where VLk 
is the angular velocity at which the star starts shedding mass through the equator. The investigation of r-modes up 
to second order in the angular velocity ft, neglecting higher-order terms arising from the expansion in a, was carried 
out in Ref. (using Saio's formalism and in Ref. [23] (using the two-potential formalism [24|). In order to neglect 
the deformation of the star due to the centrifugal force, our analysis is carried out in a slow-rotation approximation, 
ft <C VL K . 

In section HTl we review some important results obtained in the linear theory of r-modes. In section ITTT1 a nonlinear 
extension of the linear r-mode perturbation is derived up to second order in the mode amplitude a. This solution 
describes differential rotation that produces large scale drifts of fluid elements along stellar latitudes. The Lagrangian 
theory of nonrelativistic fluids |25| is used in this section to apply the boundary condition at the surface of the star and 
to show that, at second-order, r-modes do not preserve vorticity of fluid elements. The physical angular momentum 
and energy of the second-order r-mode solution are analyzed in section llVl In particular, it is shown that, in general, 
physical and canonical quantities are different. Finally, sectionals devoted to discussion and conclusions. 



II. R-MODES IN THE LINEARIZED THEORY 



The linearized fluid equations for an uniformly rotating, Newtowian, barotropic, perfect-fluid star in an inertial 
frame are 

d t S^Vi + S^v k V k Vi + v k V k 5^ Vl = -V, + <5«$) , (1) 

d t 5 (1) p + v l V l 8^ P + V l { P 5^v l ) = 0, (2) 
\7 1 V 1 6 ( - 1 ^ = 4ttGS^ p , (3) 

where v % = fWJ, and p are, respectively, the fluid velocity and the mass density of the unperturbed star, p is related 

to the pressure p by an equation of state of the form p — p(p), and the quantities S^v 1 , S^p, c^ 1 ' p, and are, 
respectively, the first-order Eulerian change in velocity, pressure, density, and gravitational potential. 

At lowest order in ft the linearized Euler equation Q yields, in spherical coordinates (r, 0, 0), the r-mode solution 

8 {1) v r = 0, (4a) 
<5 (1 V = aflCilir/R) 1 - 1 sin 1 - 1 6 sin(l(t> + cot), (4b) 
<S (1 V = u£lCil{r / R) 1 - 1 sin'" 2 9 cos 9 cos(l(f> + ut), (4c) 

with 5^U = 5 (1) p/p + <5 (1) $ given by 

<5 (1) [/ = 2aVl 2 -^—R 2 (-Y +1 sin' cos cos (/</> + wt), (5) 
I + 1 \RJ 

where u> = —Ql + 2Q/(l + 1) is the mode angular frequency in the inertial frame, R is the radius of the unperturbed 
star and Ci = (21 - 1)!! A /(2Z + 1)/[2tt(2Z)!Z(Z + 1)]. We will only consider modes for which I > 2. 

The above solution satisfies the linearized continuity equation @, which at lowest order in Q is simply V i(p5^ v l ) = 

0. 
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The linearized Euler and continuity equations do not contain information on how to split #Wf7 into the first-order 
changes in pressure S^p and gravitational potential S^Qf. This has to be done by using the perturbed Poisson 
equation J2J and the boundary conditions for the gravitational potential at the surface of the star and at infinity. 
Assuming for the same angular dependence as S^U, its (dimensionless) radial part f(r) must be a solution of 

the equation [23| 

+ im + u G /± - i±m±m m _ irG ± < n<« , (6) 

ar z r ar \ dp J dp \K/ 

with the conditions that be continuous at the surface of the star, have a continuous first derivative there and 

tend to zero as r — ^ oo. 

The r-mode solution given by eqs. JIJ and JHJ also satisfies the zero-boundary condition for the pressure p at the 
surface of the star, i.e., the first-order Lagrangian change in pressure [2^ vanishes at the surface of the star, A^p = 0. 
The Lagrangian change in pressure is, at first order, given by 

Afp = - 7 pV,£ (1) \ (7) 

where are the contravariant components of the first-order Lagrangian displacement vector and 7 = (p/p)(dp/dp). 
These components are the solution of the equations [25| 

tfdV = «9 t £ (1)i + «*v fc f (1)< - £ (1)m v m i/ = (d t + no^) (8) 
8 {1) P = -Vi{pt {x) % (9) 

where dt + 0<9^ is the time derivative in the rotating frame, dt'- At lowest order in £1 these equations yield for 
the following solution: 

£ (1)r = 0, (10a) 

= ~aCil(l + l)(^) 1 sin'" 1 6coB(l<t> + u}t), (10b) 
2 V RJ 

= \aCil(l + 1) (—) 1 sin 1 - 2 6 cos6 sm{l(t) + ujt). (10c) 

Using the above t;^ 1 it is now straightforward to verify that the Lagrangian change in pressure Ay'p vanishes 
everywhere in the star, including the surface. 

The Lagrangian displacement t;^ 1 , given by eqs. lfTU|l . is canonical |25j: it satisfies the conditions e^VjA^'vk = 0, 
where the first-order Lagrangian change in velocity is given by 



Af^ 



a t d 1) +ti*(v4 l) + v fc d 1) )- (ii) 

For the canonical Lagrangian displacement (|10|1 the canonical energy 

^^(ViVjP + pViVj®)] dV + 0{n 4 ), (12) 



1 

E < = 2 



is negative for any I > 2 and for arbitrarily slow rotation Q , implying r- modes instability to gravitational radiation |2fi| | . 
Let us emphasize here that the use of this stability criterion does not require knowledge of second-order Lagrangian 
displacements, since the canonical energy is quadratic in first-order quantities. 

Using the relations A^ ) e ijk = -e ljk V m ^ 1)m = and A^VjVk = VjA^V - v m VjV k €m, it immediately follows 
that the canonical condition e^VjAl 1 -'?;*. = is equivalent to the statement that the first-order Lagrangian change 

in vorticity is zero, A^ (e yfe VjV k ) = 0. Thus, canonical displacements are precisely those that preserve vorticity. If 
the vorticity of a r-mode perturbation is initially small, then as the perturbation grows under gravitational radiation 
reaction the vorticity will not grow. 
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III. THE SECOND-ORDER R-MODE SOLUTION 



At second order, the perturbed Euler, continuity, and Poisson equations in an inertial frame are given, respectively, 
by 

d t SWvi + S^v k V k v z + v k V k 5W Vl + 8Wv k V k 8^v t = -V^U + (^) , (13) 

d t SWp + v'V^p + V<(/xS< 2 V) + VitfMpSWv*) = 0, (14) 
V i V l <5( 2 )$ = 47rG(5( 2 V, (15) 

where S^U = S^p/p + and S^q denotes the second-order Eulerian change in a quantity q. 

Let us now assume that a S> (Q/Qk ) 2 ■ Then, terms in 8^ 2 'v l proportional to ail 3 (arising in an expansion in powers 
of the angular velocity of the star) can be neglected when compared with terms proportional to a 2 fl (arising in an 
expansion in powers of the mode amplitude). For the same reason, we neglect in S^U and S^p terms proportional 
to ail 4 . Since the second term on the right-hand side of eq. (|13fl is of order a 2 fi 4 , equation l|13|l reduces, at lowest 
order in 0, to 

d t 5 {2) v l + 8 {2) v k V k Vi + v k V k 8 {2) Vl + V 2 <S (2) [/ = -8 (1) v k V k S {1) v l , (16) 

where S^v 1 is of order aVL, S^v 1 is of order a 2 and S^U is of order a 2 Q 2 . In this equation second-order quantities 
depend on the first-order ones only through the term quadratic on S^v 1 . Using eqs. (@J to compute this term, the 
above equation reads, in spherical coordinates (r, 9, <f), 

(d t + nd <p )S {2) v r - 2VLr sin 2 6»<5 (2 V + d r 5 {2) U = 

1 / r \ 21-1 

— a 2 n 2 Cfl 2 R[-) sin 21 - 2 9 {sin 2 9 - 2 + sin 2 6»cos[2(/0 + wi)]} , (!7a) 



\R 

(d t + nd^)5 {2) v e - 2f7sin(9cos6i<5 (2 V + ^d g S (2) U = 

21-2 



^a 2 fl 2 Cfl 2 2 sin 21 - 3 9 cos 9 {sin 2 9 + 21-2 + sin 2 9 cos[2(Z</> + cut)} } , (17b) 



(d t + sw^v* + 2n— + 2n^5^v e + —L^a^u 

r sm9 r 2 sin 9 



1 / r \ 21-2 

+L a 2 n 2 c 2 i 2 ' s 

2 



(^) 2 ' 2sin2 ^ 26,sin [ 2 (^ + ^)]- ( 17c ) 



The right-hand side of the previous equations contains a piece that does not dependent on t and and a double "fre- 
quency" 2(l<fi + ojt) piece. The former induces an axisymmetric time-independent second-order solution corresponding 
to differential rotation, while the latter induces a second-order solution corresponding to an oscillating response at a 
"frequency" twice that of the r-mode. In this article we will be concerned exclusively with the solution corresponding 
to differential rotation of r-modes. 

The axisymmetric time-independent right-hand side of system (|17|) induces the following second order solution: 

5 {2) v r = 0, (18a) 
5 {2) v e = 0, (18b) 

£( 2 V = la 2 nC 2 l 2 (l 2 -I) (^-) 2 ^ 2 sin 21 -* 9, (18c) 
2 \R/ 

and 

8&U = ~a 2 n 2 C 2 lR 2 (^) 2( sin 2 '" 2 9 [sin 2 9 - 2l 2 } . (19) 

Indeed, assuming axisymmetry and time-independence system (I17|l decouples into two independent systems, one 
determining S^v^ and S^U, the other relating S^v r and S^v e . Eliminating 5^U from eqs. (|17a|l and l|17b(l . one 
obtains an equation that determines S^v^ ', namely 

fisin0cos03 r (W 2 V) - flrdg (sin 2 fl^ 2 V) = a 2 n 2 Cfl 2 (l 2 - 1)R (^) 2 ' ^ sin 2 '" 3 6>cos<9, (20) 
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yielding (|18cB . Now, inserting the above solution for S^v^ in eqs. (|17afl and (|17b|l one obtains (|19(l . Finally, through 
eq. ^7^W>v r = implies S^v e = 0. 

The above second-order solution, induced by first-order quantities, describes a drift of fluid elements along stellar 
latitudes. Note, however, that the homogeneous part of system (|17fl . involving only second-order quantities, also 
admits a solution describing a drift along stellar latitudes, 

6 {2) v r = 0, (21a) 

5 {2) v e = 0, (21b) 

S^v* = a 2 nAr N ^ sin^" 1 0, (21c) 

with 

m = 2aWA N+1 . N+le { ) 

N + 1 ' v ; 

where A and N are some constants determined by initial data. Thus, the full second-order solution describing 
differential rotation has a piece induced by first-order quantities and another one which is a pure second-order effect. 
Since the differential rotation given by eq. (|21cfl is stratified on cylinders, it cannot cancel differential rotation induced 
by first-order quantities, which is not stratified on cylinders. Thus, differential rotation is an unavoidable feature of 
nonlinear r-modes. 

The solution given by eqs. (|18|l . Ijl9|l . I|21|l and is the nonlinear extension of the linear r-mode perturbation we 
were looking for. 

The above solution also satisfies the perturbed continuity equation, which at lowest order in £1 is simply 
Vi(p6 (2 V) = 0. Again, as in the linearized theory, the splitting of S^U into 5^p and is done by using 

the perturbed Poisson equation H15[) with appropriate boundary conditions for the gravitational potential 5^$>. 

The second-order r-mode solution should also satisfy the zero-boundary condition for the pressure p at the surface of 
the star, i.e., the second-order Lagrangian change in pressure [2^ should vanish at the surface of the star, A^p = 0. 



The Lagrangian change in pressure is given by 



A 4 p 7(7 ~ 1) ( A{p\ 
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Agp = 7P— + ' P — +•••> ( 23 ) 
P 2 \ P J 

where 7 = (p/p)(dp/dp). The Lagrangian change in density A^p can be expressed in terms of the displacement vector 
using conservation of mass |25| : 

= -V fe £ (1)fe + \ (v fe e (1)fe V m £ (1)m + V fe £ (1)m V m £ (1)fe ) - V,e (2)fe + 0(e), (24) 



P 

where £( 2 ) 1 are the contravariant components of the second-order Lagrangian displacement vector. 

>? 

,,2 



Using eq. (|24|l the second-order Lagrangian change in pressure A^p is given by 



(25) 



(26) 



Since V^^ 1 ^ = 0, the boundary condition A^p — is satisfied if 
or, using eqs. (fTTHl . if 

Vfe ^ 2 ) fc = -a 2 C 2 f(l + l) 2 f -) 2 '~ 2 \{l 2 -l + l) cos 2 9 - I] sin 2 '" 4 0. (27) 
4 V R J 

The contravariant components of the second-order Lagrangian displacement vector £( 2 ) 1 are determined by the 
equations 

<5 (2 V = (dt + fld^ 2)t -£ (1)fc V fe <5 (1 V, (28) 



S^p = -p (v k ^ k - iv^^V^ 1 ^ - ^ (2)fc V fcP + ^ (1)fe C (1)m V fe V mP ) 



(29) 
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where dt + fld$ is the time derivative in the rotating frame, dt>- In the above equations, terms quadratic in the 
first-order quantities ^ 1 ' t and S^'v 1 give rise to double "frequency" terms; since in this article we are concerned only 
with differential rotation, these double frequency terms will not be considered. Thus, at lowest order in f2, eq. Il28() 
yields the solution: 

£ (2)r = d(r,0), (30a) 
C (2)e = C 2 (r,9), (30b) 
^ = D(r,9)t' + C 3 (r,9), (30c) 

where 

D{r,6) = -a 2 ntfl 2 {l + l){2l~\)(-) 2l ~ 2 s\n 2l - 2 9 + a 2 nAr N - 1 sm N - 1 9, (31) 



4 ' v '\R, 

and Ci are some arbitrary functions of r and 8. 

The second-order Lagrangian displacement vector ^'■ 2 - >I , given by eqs. I|30|) . is also a solution to eq. (|29|l at lowest 
order in fi, if the functions d are chosen to be 

1 / 7* \ 2Z— 1 

Ci(r,0) = — a 2 CfZ 2 (Z + l) 2 i?(- sin 2 '" 2 6» (sin 2 6» - 2) , (32a) 
16 vit/ 

C 2 (r,6>) = — a 2 C?l 2 (l + l) 2 f-V' 2 sin 2 '" 3 6» cos 9 (sin 2 + 21 - 2) , (32b) 
16 \RJ 

C 3 {r,9) = 0. (32c) 
It is now straightforward to verify that ^ 2 ^, given by eqs. (|30|l - (|32(l . satisfies the condition l|27[) everywhere in the 

(2) 

star, including the surface. Thus, at second-order the solution satisfies the boundary condition p = for any 
value of the constants A and N. 

At second order, the Lagrangian change in vorticity is given by 

q* = A [2 \e^ k V 3 v k ) = Afe^VjVk + e ijk Afv jVk + A^e^A^VjVk, (33) 

where 

A m e ijk = _ e ijfe Vm ^(2)m _f(i)«v„ (e iJ ' fc v m £ (1)m ) + e ljk v m (c (1)n v n ^ (1)m ) 

+e mm V m C (lb V„C (1)fc + e mj ' n V m £ (1)i V„£ (1)fc + e mnfc V m £ (1)i V n £W j (34) 

and 

Afv.-Ufc = V,-Af ^fc - v m V,-V fc ^ - (d^ (1)m + ^V^( l ' m ) VjVfc^- (35) 

Using £W* and given, respectively, by eqs. ljTU|) and eqs. (|3TJ|l - i|52|) . and taking into account that the second-order 
Lagrangian change in velocity is given by [25j 

Af Vi = d t ^ k v4 V < + v k V k ^ m V^ + d4 2) + v k (V4 2) + V,C (2) ) , (36) 
one obtains for q l the following expressions: 

q r = ia 2 fiQ 2 Z 2 (Z + l)(-^) 2Z ~ 2 sin 2i - 4 6lcos6'[;(/ + l)sin 2 6' + 4(;-l) 2 ] 

+a 2 ilA(N+ l)r N - 1 sm N ~ 1 9cos9, (37a) 
--a 2 nCfl 3 {l + l)iT 1 (-) 2l ~ 3 sin 2 '- 3 9[(l + 1) sin 2 9 + 4(1- 1)] 



q = 4 __._ 1 ^ 

-a 2 QA{N+ l)r N ~ 2 sm N 9, (37b) 
g* = 0. (37c) 

As it can be seen from eqs. (|37|) . the second-order Lagrangian change in vorticity is always different from zero for 
any values of the constants A and N. Thus, at second-order, r-modes do not preserve vorticity of fluid elements. 

As it was mentioned in section [D] the first-order canonical displacements £W* are precisely those that preserve 
vorticity. As a result, if the vorticity of a linear r-mode perturbation is initially small, then as the perturbation grows 
under gravitational radiation reaction the vorticity stays small. However, the situation changes when one takes into 
account second-order effects; due to differential rotation an initially small vorticity may increase as the perturbation 
grows. 
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IV. ANGULAR MOMENTUM AND ENERGY OF THE R-MODE PERTURBATION 



The physical angular momentum of the second-order r-mode solution found in the previous section is, at lowest 
order in f2, given by 

5 {2) J = J pSWvfdV 

= -(l-l)(2l + l)a 2 rtR 2 - 21 [ R pr 2l+2 dr 
2 Jo 

+ 2ira 2 nA / pr N+3 dr / sin N+2 6d6, (38) 

Jo Jo 

where S^v^ is given by eqs. I|18cfl and i|21c|l . 

Since the physical angular momentum depends on the arbitrary constants A and N, one would like, at this point, 
to specify some physical conditions that fix the values of these constants, i.e., that pick out a member of the family 
of the second-order solution given by equations (|18|l . I|19|) . (121(1 . and Since in the linearized theory vorticity is 

conserved, it would be natural to impose the condition that the Lagrangian change in vorticity at every order should 
be zero. However, as we have seen, due to the presence of differential rotation, vorticity is not conserved already at 
second order. 

The second-order physical angular momentum can be decomposed in two pieces [25l | ; one linear in the second-order 

(2) 

Lagrangian change in velocity Vi and another one, called the canonical angular momentum J c , quadratic in the 
first-order Lagrangian displacement vector 

5(2)j = 77 / pv l A {2) Vl dV + J c , (39) 



where 

J, = - 



P d^ 1 [d t tf ] + « fc v fc & {1) ) dv (40) 



and the second-order Lagrangian change in velocity is given by eq. (|36[) . 

For the Lagrangian displacement given by eqs. I|10l) and eqs. (|30[) -H32 [I . the canonical angular momentum and the 
integral of the second-order Lagrangian change in velocity are computed to be, respectively, 

rR 



I r 
J c = --l{l + l)a 2 nR 2 - 21 \ 
4 Jo 



pr M+2 dr (41) 



and 



i J pv l A {2) Vl dV = i (hi 2 -1-2) a 2 nR 2 - 21 j\r 2l+2 dr 

+ 2ira 2 nA / pr N+3 dr / sm N+2 6d6. (42) 
Jo Jo 

As it can be seen from the above equations, the second-order physical change in the angular momentum is not 
equal, in general, to the canonical angular momentum. Since the second-order physical change in the energy (in the 
inertial frame) S^ 2 'E and the canonical energy E c are related by the expression [25j 

6WE= { pv l /S. [2) v,dV + E C1 (43) 



5( 2 'E and E c are also, in general, not equal. 

In recent years, the nonlinear behaviour of r-modes was investigated within a toy model of a spherical shell of rotating 
incompressible fluid |20| . It was shown that in a spherical shell r-modes carry zero physical angular momentum, 
£( 2 ) J = 30], and positive physical energy, S^E > 0, while both the canonical angular momentum and canonical 
energy are negative. Thus, for r-modes in a spherical shell one cannot equate physical and canonical quantities. Based 
on these results, it was conjectured in Ref. [2(J that also in the case of a full star physical angular momentum and 
energy of r-modes are not equal to the corresponding canonical quantities. Our investigation, on r-modes at second 
order in a, confirms that for a full star, in general, 5&> J ^ J c and ^= E c . Note, however, that for specific choices 

of initial data [such that the two terms on the right hand side of eq. (|42|) cancel each other] physical and canonical 
quantities can be made equal. 
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V. DISCUSSION AND CONCLUSIONS 

We have found a nonlinear extension of the linear r-mode perturbation, describing differential rotation of pure 
kinematic nature that produces large scale drifts along stellar latitudes. Differential rotation of fluid elements given 
by eqs. (|18|l is induced by first-order quantities, while differential rotation given by eqs. i|21[l is a pure second-order 
effect. The latter is stratified on cylinders and, therefore, cannot cancel differential rotation induced by first-order 
quantities, which is not stratified on cylinders. As already mentioned, our computation was carried out in a slow 
rotation approximation, £1 <C Qk', in order to neglect higher-order terms arising from the expansion in powers of 
the star's angular velocity f2, it was also assumed that a ^> (VI/VIk) 2 [see discussion before eq. Recent results 

indicate that in rapidly rotating neutron stars (just born in a supernova or spun up by accretion in low-mass X-ray 
binaries) the r-mode instability may saturate at low values [lfl Il7| , implying that a could be of the same order of 
magnitude as (il/flfc ) 2 . In that case, terms in 8^v l proportional to aQ 3 (arising in an expansion in powers of the 
angular velocity of the star) are as important as terms proportional to a 2 fl (arising in an expansion in powers of 
the mode amplitude) and, therefore, should be taken into account. However, even though the velocity field 5^v % 
and other second-order quantities derived in this paper are, strictly speaking, only valid in the regime f2 <C and 
a ^ (il/ilx) 2 , they could be used to illustrate the influence of differential rotation of r-modcs on the nonlinear 
evolution of rapidly rotating neutron stars. 

Recently, an analytical expression for the azimuthal drift velocity of the / — 2 mode was derived from the linearized 
fluid equations by expanding the velocity of a fluid element located at a certain point in powers of a, averaging over 
a gyration, and retaining only the lowest-order nonvanishing term |l7j |. A comparison with our results shows the 
expression obtained in Ref. [l7j is qualitatively correct but not exact to 0(a 2 ). This is to be expected since the 
procedure used there does not consider nonlinear effects in the fluid equations. 

In the linearized theory, r-mode perturbation preserves vorticity of fluid elements, with the consequence that the 
vorticity will not grow as the perturbation grows under gravitational radiation reaction. However, the situation 
changes when one takes into account second-order effects; due to differential rotation of fluid elements, producing 
large scale azimuthal drifts, the second-order Lagrangian change in vorticity is always different from zero and may 
increase as the perturbation grows under gravitational radiation reaction. 

It was also explicitly shown that for r-modes physical and canonical quantities cannot be equated. Canonical angular 
momentum (or energy) is not the full angular momentum (or energy) at second order; one should also include a part 
linear in the second-order Lagrangian change in velocity, which, as pointed out in Ref. pfij . is related to conservation 
of circulation in the fluid. Since, as shown above, at second-order r-modes do not conserve vorticity, it follows that, 
in general, the physical and canonical angular momentum (or energy) do not coincide. However, specific choices of A 
and N can be made such that the integral in the right-hand side of eqs. I|39fl and (|43|l vanishes. Such a case (for I = 2) 
was studied in Ref. within a phenomenological model for the evolution of r-modes. However, it is not clear which 
physical condition would force the constants A and N to take such particular values; thus, it seems more appropriate 
to study the evolution of r-modes for the case of arbitrary values of A and N. 

In recent years, numerical simulations have been used to study the nonlinear evolution of r-modes, both for rela- 
tivistic and Newtonian stars. Initial data for this evolution was generated by adding a linear perturbation 5^v l to 
an equilibrium stellar model. An analytical expression for <5W v l is known only for small-amplitude perturbations of 
slowly-rotating Newtonian stars. As shown in this paper, the linear r-mode perturbations given by eq. JIJ induce, 
at second order in the mode amplitude, a drift of fluid elements along stellar latitudes, which cannot be avoided by 
special choices of the constants A and N. Thus, initial data for numerical non-linear evolution of r-modes should also 
include a piece describing differential rotation. 
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